Correlations of electromagnetic fields in chaotic cavities 
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We consider the fluctuations of electromagnetic fields in chaotic microwave cavities. We calcu- 
late the transversal and longitudinal correlation function based on a random wave assumption and 
compare the predictions with measurements on two- and three-dimensional microwave cavities. 



Classical ergodicity suggests that wave functions in 
chaotic systems may be described by superpositions of 
waves with wave vectors of constant length but random 
directions E|,0|. Their fluctuations are distinctly differ- 
ent from the more familiar optical speckle patterns where 
also the wave numbers fluctuate || . The distributions of 
amplitudes turns out to be Gaussian and the spatial auto- 
correlation function is given by Bessel functions of order 
| — 1, where d is the billiard dimension. Overwhelm- 
ing evidence for this has been accumulated especially in 
numerical studies of billiards [^-^), experiments on mi- 
crowave billiards and surfaces of constant negative 
curvature Higher order moments and correlations 
pO| , p| as well as contributions from prominent classi- 
cal structures ('scars') have also been studied |1S^ ], so 
that by now one has a fairly good understanding of the 
fluctuation properties of scalar wave functions in chaotic 
systems. 

Recent experiments on the elastodynamics of vibrat- 
ing blocks |l3j and 3-d microwave resonators |lj,|l5| deal 
with situations where the wave fields have more than 
one component which are typically mixed by the bound- 
ary conditions |l6|-|l9|| . New effects such as ray splitting 
p8| or chaotic features in systems with integrable ray 
dynamics [fl9| are found. 

We are here interested in the consequences of these 
additional degrees of freedom for the fluctuations of the 
wave functions, fn particular, we will show that the fact 
that in the absence of space charges electromagnetic fields 
are divergence free implies differences between the longi- 
tudinal and transversal correlation functions and devia- 
tions from the behaviour expected for scalar fields || . We 
present results for the field components, the intensities 
and the frequency shift and compare with experiments 
on microwave billiards. 

Our starting point is the electromagnetic analog of the 
the semiclassical ansatz that the scalar field is a superpo- 
sition of plane waves with constant amplitudes and fixed 
wave length but randomly oriented wave vector For 
the electromagnetic case we have in addition to allow for 
different orientations of the polarization. We thus as- 
sume that the field at a point x in position space is due 
to a superposition of many plane waves with uniformly 
distributed orientations of polarization and wave vector, 
e.g., 

E(x)«^E„ e * k - x , (f) 
and similarly for the B-field. The complex amplitudes 



~E U and of all waves are transversal, E„ ■ k„ = and 
B^ • k u — 0, and satisfy E„ • B y = so that k„, E„ 
and B, y form an orthogonal dreibein. The absolute val- 
ues |E„| and |B, y | are all the same, that is to say, we 
assume that there are no losses during reflections at the 
walls. From this assumption there follows immediately 
that all three components are Gaussian distributed with 
the same distribution p"4[ ] 

The spatially averaged correlation functions then are 
for the electric field 

C B ,y(r) = (E t (x + r/2)E,(x-r/2))/(E 4 2 ), (2) 

and similar for the magnetic field and the cross correla- 
tion between E and B. The normalization is by the mean 
square amplitude of a single component of the fields, 
(Ef), so that Ce,u{0) = 1- In a tensor notation, this 
can be combined to a tensor of correlation functions, 

C E (v) = (E(x + r/2) ® E(x - r/2)}/(E?) . (3) 

Substituting ([!]) and performing the spatial average then 
results in (the normalization will be restored in the end) 



C E (r) a^E^E^ 



(4) 



To proceed further, let r point in the z-direction and 
introduce spherical coordinates for the wave vector, 



(5) 



The electromagnetic field contributions lie in a plane per- 
pendicular to this wave vector, spanned by the two vec- 
tors 
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If i/v denotes the angle of polarization, the field compo- 
nents are 



E„ = cos i\) v + sin ip v 
By = — sini/v + cos if) v 



(7) 
(8) 



In the limit of a large number of contributing compo- 
nents, the sum over the different contributions can be 



1 



replaced by a continuous average over all directions (an- 
gles 9 and <fi) for the wave vector and all polarizations 
(angle V), 
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After averaging over the polarizations and the azimuthal 
angle, the correlation functions become 
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C E (r) = - 




ikr ( 



(10) 



The final average over 9 can be expressed in terms of 
spherical Bessel functions, but it is more convenient to 
use the trigonometric representation directly, 




//x(fcr) 
C E (r) = /x(fcr) 
V 



with the transversal correlation function 



3 / sin£ sin£ - £cos£ 
/±(0 = 2 ? 



and the longitudinal correlation function 
sin £ — £ cos £ 



(11) 



(12) 



/||(0 = 3- 



The asymptotic behaviour of these functions is that for 
small r they approach the same value, Ce,u — > 1, by nor- 
malization. For large r they oscillate on a scale set by 
the wavenumber and decay like 1/r for the transversal 
and like 1/r 2 for the longitudinal correlations. For the 
trace of the correlation function, 



trC f 



(E(x + r/2) -E(x- 
_ sin kr 



r/2))/(E 2 



kr 



(14) 
(15) 



the correlation between polarizations and wave vector is 
eliminated and one recovers Berry's result for random 
waves in three dimensions || , except for a factor due to 
the normalization (see Fig. 1). 



FIG. 1. Correlation functions for longitudinal and 
transversal fluctuations (dashed and dotted line, respectively) 
for the electromagnetic fields. The full line shows the correla- 
tion function for a 3-d scalar wave function for comparision. 

The correlations for the magnetic field B have the same 
functional dependence, 



C B (v) = 




(16) 



There are no correlations between E and B. 

For the experiments also the correlations of intensities, 



C EE (r) = (|E(x + r/2)| 2 |E(x - r/2)| 2 )/(|E| 4 



>| 4 ) (17) 



and similarly for the magnetic field are relevant. For ease 
of comparison with numerical data, we normalize this 
function by the second moment of the intensity so that 
C — > 1 as r — » 0. For large r the correlations between 
intensities decay and the correlation function approaches 
(13) (|E| 2 ) 2 /(|E| 4 ). For Gaussian random field components, 
this ratio is 3/5. Therefore, the correlation function with 
the above normalization becomes 



Cee{v) = 1 [h(kr)f + -| [/||(fcr)] 2 + ^ (18) 



and similar for the magnetic field. For the intensity cor- 
relation function of a 3-d scalar field the correlation func- 
tion becomes 
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3' 



(19) 



where the asymptotic value of 1/3 reflects the ratio of 
the square of the second moment to the fourth moment 
as expected for a Gaussian distribution MfyJ. 
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FIG. 2. Comparison between experimental results for spa- 
tial autocorrelation functions in a quarter stadium microwave 
billiards and theory. The experimental curves were obtained 
by superimposing the results for the 20 lowest lying eigenres- 
onances. (a) Spatial autocorrelation function of the electric 
field amplitudes (full line) and theoretical predictions from 
Eq. (poj). (b) Spatial autocorrelation function for the squares 
of the electric field amplitudes using the same data set as 
in (a). The dashed line corresponds to the prediction from 
Eq. @. 

To test these predictions we measured the field dis- 
tributions in two- and three-dimensional microwave bil- 
liards. We start with the discussion of the results in a 
resonator of the shape of a quarter stadium billiard. We 
measured the microwave transition amplitudes between 
two antennas, one kept fixed, the other moved around to 
probe the spatial distribution of the wave functions. At 
an eigenfrequency such measurements yield directly the 
electric field strength E(x) as a function of the position. 
Details of the experiment are described elsewhere [ pp[ . 
For microwave frequencies below v max = c/2d, where d 
is the height of the resonator, only TM modes are excited 
and the electric field has a single component E z . For this 
component, Berry's arguments give a spatial autocorre- 
lation function 

C E (r) = (E,(x + r/2) • E z (x - r/2)) ~ J (fcr) (20) 

Fig. 2(a) shows the experimental autocorrelation func- 
tion, normalized to Cb(0) = 1. It was obtained by su- 
perimposing the results from the 20 lowest eigenfrequen- 
cies of the quarter stadium. Apart from a discrepency 
of about 10 percent in the wavelength of the oscillations 
the experiment reproduces the prediction by Eq. ( po|) per- 
fectly. 

The autocorrelation function of the field intensities be- 
comes in the 2-d case 

C EE (r) = (\E z (x + r/2)| 2 \E z (x - r/2)| 2 ) 



~llMkr)f + l. (21) 

This correlation function, which has also been studied 
by Sridhar et al. Q, is shown in Fig. 2(b) for the same 
data set that entered Fig. 2(a). As in Fig. 2(a) we note 
a small difference in the wavelength of the oscillations. 
This has not been observed in the experiments of Srid- 
har et al. Q and in the stadium wave functions in Ref. [^|, 
but it has appeared for wave functions of an octogon bi - 
liard on a surface with constant negative curvature [^|, 
where it has been attributed to anisotropics in the wave 
functions. In the present case the discrepancy may be 
caused by higher order corrections to the semiclassical 
predictions since the wave functions are not very far into 
the semiclassical regime: Typical wave lengths for the 
wave functions that entered the analysis of Fig. 2 are 
about 0.2-0.5 stadium widths. Assuming corrections to 
be of order 1/fc 2 with a prefactor of order 1, as suggested 
by the analysis of plfl , the deviations can be estimated 
to be about 10%, as observed. 
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FIG. 3. Spatial autocorrelation function CAv{kr) for the 
frequency shifts obtained by the perturbing bead method in a 
three-dimensional Sinai billiard. The dashed line corresponds 
to the theoretical prediction from Eq. ^^), taking into ac- 
count the vector properties of the electromagnetic fields, (a) 
Correlation function on a linear scale, (b) The correlation 
function C& v {kr) — C& v (oo) on a logarithmic plot. The dot- 
ted line corresponds to Eq. ( |l9| ) and would have been expected 
for three-dimensional scalar fields. 

Turning to the three-dimensional microwave cavities 
note that Maxwell's equations can no longer be reduced 
to a scalar wave equation, so that effects due to the vector 
properties of the electromagnetic field can become essen- 
tial. The field distributions in a cavity of the shape of a 
three-dimensional Sinai billiard were mapped by means 
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of the perturbing bead method [Q . The technique uses 
the fact that a spherical metallic bead in the resonator 
shifts the eigenfrequency of a resonance by an amount Aj/ 
proportional to — 2E 2 + B 2 , where E and B are the fields 
at the position of the bead. Since the electric and the 
magnetic field components are uncorrelated, the spatial 
autocorrelation function for the frequency shift, 
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C Au (r) = (Ai/(x + r/2) • Ai/(x - r/2)} 



(22) 



is given by the intensity autocorrelation function de- 
fined in Eq. (|1S|), up to an off-set resulting from the 
fact that C All (r) does not vanish in the limit r — > oo. 
This off-set is easily calculated. Using again that for 
Gaussian distributions the fourth moment amounts to 
three times the square of the second moment, we find 
Cai/(0)/Cai/(oo) = 13/3. After normalization we thus 
have 

C Av {r) = § IfAkr)} 2 + ^ [/||(fcr)] 2 + ^ (23) 

The experimental results in Fig. 3(a) are in good agree- 
ment with the theoretical prediction from Eq. (^3|) for 
Ax-values below 6. The hole in the histogram for kr < 0.2 
reflects the minimum grid size used in the bead measure- 
ment. In Fig. 3(b) the data for C& v (r) ~ C Aiy (oo) are 
replotted on a logarithmic scale to emphasize the min- 
ima. Also shown is the scalar correlator (19). Note that 
the scalar function has zeroes but the spectral correlator 
does not. The absence of zeroes or, in the case of finite 
resolution, deep minima, in the experimental data is thus 
clear evidence for the influence of the polarizations. 

In summary, we have shown that longitudinal and 
transversal correlation functions of electromagnetic 
waves in chaotic cavities differ because the waves are 
transversal. The experiments verified the effect for the in- 
tensities, and more direct studies of the fields themselves 
are desirable. We expect that also in other situations 
with non-scalar waves fields, as for instance in acoustics 
in anisotropic media or hydrodynamic waves, the corre- 
lations will be characterized by tensors which depend on 
the character of the modes and the directions. 
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